EIGENVALUE ASYMPTOTIC OF ROBIN LAPLACE OPERATORS ON 
TWO-DIMENSIONAL DOMAINS WITH CUSPS 



HYNEK KOVARIK 



Abstract. We consider Robin Laplace operators on a class of two-dimensional domains 
with cusps. Our main results include the formula for the asymptotic distribution of the 
eigenvalues of such operators. In particular, we show how the eigenvalue asymptotic de- 
pends on the geometry of the cusp and on the boundary conditions. 
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1. Introduction 

Let O C be an open domain such that the spectrum of the Dirichlet Laplacian — 
on Q is discrete. Denote by Nx{—Aq) the counting function of — A^, i.e. the number of 
eigenvalues of — A^ less than A. The classical result by H. Weyl, [We| . states that if $7 is 
bounded, then 

Nx{-AE) = ^M+oi\) A^oo, (1.1) 

where \^}\ denotes the volume of il.. The proof of for unbounded domains with finite 
volume is due to M. Birman, M. Solomyak and B. Boyarski, see e.g. [BiSoj . The situation 
is different for the Neumann Laplacian — A^ . In this case equation (II. ip . with N\{—Aq) 
in place of Nx{—A^), holds whenever O is bounded and has sufficiently regular boundary, 
see e.g. [IvH [Nl INSj for the estimates on the rest term in (jl.ip . However, the Neumann 
Laplacian might not satisfy (II. ip (its spectrum might even not be discrete) if 17 has rough 
boundary or if Q is unbounded, [Ml EH iHSSl iJMSl iNSllSoI] . 

Here we will focus on unbounded domains with regular boundary and we will consider 
two-dimensional domains of the form 

n = {ix,y) GR^ : x>l,\y\< fix)}, (1.2) 

where / : (1, oo) — )■ M is a positive function such that f{x) — )■ as j; — )■ oo. Then the counting 
function Nx{—A^) of the Dirichlet Laplacian satisfies (11. ip as long as / is integrable. If / 
decays too slowly, so that \^\ = oo, then the spectrum of — A^ is still discrete, but Nx{—Aq) 
grows super-linearly in A, see \Be\ \Da\ \Ko\ ISij . On the other hand, the spectrum of the 



Neumann Laplacian — Aq is discrete if and only if 



hm ( / ^1 ( / f{t)dt]=0. (1.3) 



x—^oo 



fit) 



Date: June 30, 2010. 



1 



2 



HYNEK KOVARIK 



This remarkable fact was proved in |EHj . see also [Maj . Asymptotic behaviour of Nx{—Aq) 
on domains of this type was studied in |Bel EH [JMSl El [Sol]. We would like to point out 
that / must decay faster than any power function for (jl.3p to hold. We thus notice a huge 
difference between the spectral properties of — Aq and — Aq on such domains. 

Motivated by this discrepancy, we want to study the gap between Dirichlet and Neu- 
mann Laplacians. To do so we consider a family of Laplace operators on J7 which formally 
correspond to the so-called Robin boundary conditions 
du 

— {x,y) + h{x)u{x,y)=0, x>l, y = ±f{x), (1.4) 

where ^ denotes the normal derivative of u and h : (l,oo) — )• M-)_ is a sufficiently smooth 
bounded function. The extreme cases h = and h = oo correspond to Neumann and 
Dirichlet Laplacians respectively. First question that arises is under what conditions on h 
and / is the spectrum of the associated Robin Laplacian discrete. Next we would like to 
know how the coefficient h{x) of the boundary conditions affects the asymptotic distribution 
of eigenvalues of the Robin Laplacian. 

The paper is organised as follows. In section[3]we formulate our main results, see Theorems 
13.31 and 13.61 Similarly as in [Berl IJMSj , we show that the leading term of the eigenvalue 
asymptotic has two contributions, one of which results from an auxiliary one-dimensional 
Schrodinger operator. The boundary conditions affect the eigenvalue asymptotic through the 
term h{x)y^l + f'[xY/f{x) which enters into the potential of this operator, see equations 
(|2.2p and (|2.7|) . For some particular choices of h and / this contribution can be calculated 
explicitly, the corresponding results are given in section 13.11 

The proofs of the main results are given in section [3 Our strategy is to treat separately 
the contribution to the counting function from a finite part of Vl and from the tail. In 
section 15.11 it is shown that the contribution from the finite part satisfies the Weyl law 
(jl.ip . The key point of the proof is to transform, in the remaining part of the problem 
to a Neumann Laplacian plus a positive potential that reflects the boundary term, see 
section 15.21 To this end we employ the technique known as ground state representation, 
which has been recently used e.g. in [FSWj to derive eigenvalue estimates for Schrodinger 
operators with regular ground states, see also [FLSj . Once this transformation is done, we 
show, by rather standard arguments, that one part of the eigenvalue distribution of such 
Neumann Laplacian with additional potential is asymptotically (i.e. for A — ?• oo) equivalent 
to eigenvalue distribution of a direct sum of certain one-dimensional Schrodinger operators, 
see section 15.31 This enables us to prove Theorem 13.31 Finally, in the closing section [6] we 
discuss some generalisations for Robin Laplacians with non symmetric boundary conditions. 

2. Preliminaries and notation 

Given a self-adjoint operator T with a purely discrete spectrum we denote by Nx{T) the 
number of its eigenvalues, counted with multiplicities, less than A. We will write A c^l B \i 
the operators A and B are unitarily equivalent and we will use the notation 

/i(A)~/2(A) A^oo ^ iim:M^ = l. 

A^-oo j2(Aj 

We will consider the eigenvalue behaviour of the Robin boundary value problem in a weak 
sense. Therefore the main object of our interest is the self-adjoint operator A^j in L^(J7) 
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associated with the closure of the quadratic form 

Q„[u]= [ \Vu\'^dxdy+ [ a{x){\u{x,f{x))\'^ + \u{x,-f{x))\'^)dx (2.1) 
Jn Ji 

on Cq(S7). Here Cq(0) denotes the restriction to of functions from C^(M^) such that 
for each y the support of u{-,y) is a compact subset of (l,oo). The operator A^^ formally 
corresponds to the Laplace operator on Q with Dirichlet boundary condition at {x = 1} and 
mixed boundary conditions (II. 4p at the rest of the boundary, if we chose a such that 

a{x) = h{x) y^lTTW . 

Remark 2.1. Since we work under the assumption that f'{x) — )■ as x — )■ oo, see below, 
and since the asymptotic of Nx{A„) depends only on the behaviour of a at infinity, from 
now on we will work with the function a instead of h. 

We will also need the following auxiliary potentials: 
Throughout the whole paper we will suppose that / satisfies 

Assumption 2.2. / G C°°(l,oo) is positive and such that f'{x) < for all x large enough. 
Moreover, 

lim f{x) = lim f"{x) = 0. (2.3) 
Note that ()2.3p implies /'(x) — )• as x — )• oo, see Lemma [23] below. 

Lemma 2.3. Let f G C^(l,oo) be a nonnegative function. Assume that f and \ f"\ are 
bounded on (l,oo). For a given x > 1 define = sup^>^ /(s) and M'J. = sup^>^ |/"(s)|. 
Then 

if'{x)f<2M,M'^. (2.4) 
Proof. Let s > x. The Taylor expansion of / at the point x gives 

/(s)-/(x) = t/'(x) + ^r(y), yG[x,s], t = s - x. (2.5) 

On the other hand, / > ensures that |/(x) — f{s)\ < for all s > x. This together with 
(|2.5p implies that the inequality 

holds for all t > 0. Optimization with respect to t then gives the result. □ 

Remark 2.4. Note that if we leave out the assumption / > 0, then the above proof still 
works with the modification that now |/(x) — f{s)\ < 2Mx. This results into the Landau 
inequality, i.e. inequality ()2.4p with the factor 2 on the right hand side replaced by 4. 

The hypothesis on a are the following: 

Assumption 2.5. The function a G C'^{1, oo) is non negative. Moreover, a, a' and a" are 
bounded and 

lim W„{x) = oo. (2.6) 

a;— ^oo 
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In order to formulate our next assumption, we introduce the operator 

'Ha = --r^ +W^{x) in L2(1,oo) (2.7) 

with Dirichlet boundary condition at x = 1. More precisely, is the operator generated 
by the closure of the quadratic form 

(iV'f + V'^) dx, ^ G Cl{l, oo). 
Alongside with T-L^j we will also consider the auxiliary operator 

^ = -^--J^x)^y i>'((l,oo)x(-l,l)) (2.8) 

with Dirichlet boundary conditions. 
Assumption 2.6. For < e < 1 we have 

Nx{{l±e)n^) = Nx{n^){l + 0{e)), (2.9) 
N^{{l±e)B) = N^{B){l + 0{e)) (2.10) 

Remark 2.7. A similar assumption was made in |JMSj . Although this assumption is es- 
sential for the approach used in the proof of Theorem 13.31 below, it is natural to believe that 
the statement holds under more general conditions. Note also that for domains with finite 
volume (|2.10p holds automatically. 

3. Main results 

Theorem 3.1. If \2.S\ and \2.5\ are satisfied, then the spectrum of A„ is discrete. 

Remark 3.2. Contrary to the case of Neumann Laplacian, the spectrum of A„ can be 
discrete also if the volume of is infinite. For example if a is constant, then (|2.6p is 
automatically satisfied in view of the fact that f{x)V{x) — )■ as x — )■ oo, see equation (j5.13p . 
On the other hand, condition (12. 6p is, unlike (II. 3p . only sufficient. 

Theorem 3.3. Suppose that assumptions \2.^ \2.5\ and \2.6\ are satisfied. Then 

Nx{A^) ~ Nx{-AE) + Nxin^) A ^ oo. (3.1) 

Remark 3.4. The second term in ()3.ip is a contribution from the eigenvalues of the operator 
A„ restricted to the space of functions which depend only on x. This is analogous to the case 
of Neumann Laplacian, |DS1 IJMSl ISol] . On the other hand, the presence of the boundary 
term cr{x) enables us to apply ()3.ip also in the situation in which the Neumann Laplacian 
does not have purely discrete spectrum. 

Remark 3.5. Theorem 13 . 3 1 allows a straightforward generalisation to Robin Laplacians with 
different boundary conditions on the upper and lower boundary of il, say given through 
functions (Ti(x) and (72{x). In that case we only have to replace a{x) in (|2.2p by (o"i(x) + 
(^2{x))/2, see section [6TT] for details. 

For domains with finite volume Theorem 13.31 and the Weyl formula (jl.ip give 



EIGENVALUE ASYMPTOTIC OF ROBIN LAPLACE OPERATORS 



5 



Theorem 3.6. Let \Q\ < oo and suppose that assumptions \2. Si \2.5\ and ()2.9p are satisfied. 
Then 

Nx{A^) ~ ^ l^^l + Nxi-Ha) A ^ oo. (3.2) 
47r 

Remark 3.7. Note that if a = 0, then the condition < oo is necessary for the spectrum 
of Af] = — Ao to be discrete, see (11.31). Hence in that case there is no difference between 



Theorems 13.31 and 13.61 and the resulting eigenvalue asymptotic agrees with the one obtained 
in [B^lJMS] . 

Corollary 3.8. Let \Q\ < oo and let o'{x) = a be constant. Assume that f satisfies \2.SX 
Then 

limsup 2;V(a:) = =^ Nx{A„) \ A ^ oo (3.3) 



\n[ \a\ \ 



lim x'fix) =a' =^ Nx{A^) ^ { ^ + ] X A^oo. (3.4) 

x^oo \An 4yo" / 

Remark 3.9. Equation (13. 3p provides a sufficient condition on the decay of / for the Weyl's 
law in the case of constant a. Notice that the borderline decay behaviour is /(x) ~ 
which is in contrast to f{x) ~ x~^ in the case of Dirichlet Laplacian. The reason behind 
this is that the principle eigenvalues of Robin and Dirichlet Laplacians on an interval of the 
width 2f{x) scale in a different way as f{x) — )• 0. Observe also that (|3.4p turns into (|3.3p 
when fj — )• oo, as expected. 

If the volume of is infinite, then we confine ourselves to situations when / is a power 
function. The asymptotic distribution of the Dirichlet-Laplacian on such region is known, 
see |Roj , [Sij . These results together with Theorem 13.31 yield 

Corollary 3.10. Let f[x) = x~°', < a < 1. LfHa satisfies (|2.9p . then as A — )■ oo we have 
vr 

where C(-) and B{- , ■) denote the Riemann zeta and the Euler beta function respectively. 
3.1. Examples. We give the asymptotic of N\{A„) for some concrete choices of / and a. 
3.1.1. f{x) = a > 1, (t(x) = a = const. Here 

Wa{x) = ( Y + f ) + 

is convex and increasing at infinity so that assumption 12.61 is satisfied, see |TH Chap. 7]. 
Theorem 13.61 in combination with Theorem 14. 2| see Section HI gives 

ArA(^a) ~ T^A + — ( i, ^ ) A5+i A^oo. (3.5) 

An an \a z J 

Note that, in agreement with Corollarv 13.81 -^a(^o-) obeys Weyl's law as long as a > 2 and 
for a = 2 the order of A is linear, but the coefficient is different from the one in the Weyl 
asymptotic. When a < 2, then the behaviour of N\[A^) for A — )• oo is fully determined by 
the second term on the right hand side of p.Sp . 
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3.1.2. f{x) = < ce < 1, a{x) = ax~^. Assumptions 12.51 is satisfied if and only if 
< 13 < a. For these values of /3 Corollary 13.101 and Theorem 14.21 give 

1 

(a — pjTT \a — p 2 J 

4. Auxiliary material 

In this section we collect some auxiliary material, which will be used in the proof of the main 
results. First we fix some necessary notation. Given a continuous function q : (l,oo) — >• M 
such that q{x) — )• oo as x — )• oo, we denote by T^^'^.^ the operator in L^(o, b) acting as 

= + 1 < a < 6 < oo 

with Dirichlet boundary conditions at x = o and x = b. Operators T^^'^ , T^'^ and T^'l^ 
are defined accordingly. For 6 = oo we use the simplified notation etc. to indicate 

the corresponding boundary condition at x = a. It is well known that imposing Dirichlet 
boundary condition at x = a is a rank one perturbation. Variational principle thus implies 
that 

0<iVA(T£oo))-^^A(r£oo))<l Va. (4.1) 

Lemma 4.1. Suppose that q{x) is a continuous function such that q[x) — ?• oo as x ^ oo. 
Then for any s > 1 it holds 

MTiloo)) ~ MTiloo)) ~ MtS,oo)) - M^iloo)) A ^ oo. (4.2) 

Proof. In view of (j4.ip it suffices to consider the Dirichlet operator only. Let Ix ■= {x > s : 
q{x) < A/2}. Then 

where the superscript Dir indicates Dirichlet boundary conditions at the end points of Ix. 
Since — )• oo as A — ?• oo, this shows that liminfAH>oo ^^"^^"^ -^>'i'^l^ oo)) — view of 

the equation 

Mt;I:^)) - MT^If)) = 0{VX) A ^ oo V. > 1, 
the result follows from the Dirichlet-Neumann bracketing (by putting additional boundary 
conditions at x = s), see e.g. \RS\ Chap. 13]. □ 

Under certain additional assumptions one can recover the eigenvalue distribution of such 
operators from the potential q. The following theorems are due to [Til Chap. 7]: 

Theorem 4.2 (Titchmarsh). Suppose thatq{x) is continuous increasing unbounded function, 
that q'{x) is continuous and x^q'{x) — )■ oo as x —)• oo. Then 

(t£oo)) ~ - y (A - qixM dx, A ^ oo. (4.3) 

Theorem 4.3 (Titchmarsh). Suppose that q(x) is continuous increasing and convex at in- 
finity. Then (j4.3p holds true. 

A simple combination of the above results gives 
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Lemma 4.4. Assume that there exists some Xc such that q satisfies the hypothesis of The- 
orem \4^ or \4 -31 for all x > Xc- Then for any s > Xc we have 

1 />oo 1 1 poo I 

Nx (T(?^)) ~ - (A - q{x))l dx- - (A - q{x))l dx A ^ oo. (4.4) 

Next we consider the operators 

B^/^ = -dl - d'^ in L2((n,oo)x(-l,l)) 

subject to Dirichlet boundary conditions on (n, oo) x ({1} U {—1}) and Neumann/Dirichlet 
boundary condition on {n} x (—1,1) respectively. We have 

Lemma 4.5. For any n G N holds 

A f'°° 

N^{B^) ^ N^{B^) ^ — j f{x)dx A^oo if < oo, (4.5) 

{B^) ~ iVA (Bf ) ~ N),{B) A ^ oo if M= oo. (4.6) 

Proof. Equation (j4.5p for B^ fohows directly from }SVi Thm. 1.2.1]. Hence it remains to 
prove ^M) for B^ and (gj]). Note that 



d^ Tx^e 



N,{B^) = Y,mL^^J, iV,(^^) = J]iV,(L^J, where L^/f = - + 

k=l k=l ' 

are one-dimensional operators acting in L'^{n, oo) with Neumann/Dirichlet boundary condi- 
tions on X = n. Obviously there exists a positive constant c such that for any n and any k 
the operator inequality L^^ > > ck"^ holds. This means that there exists some K{X) 
with K{X) = 0(\/A) as a — )• oo and such that 

A-(A) KiX) 

j;iV,(L^J = iVA(L^J, E^a(L^J = E ^A(^^n) 

fc>l fc>l fc>l fe>l 

Moreover, since < Nx{L^^J - NxiL^J < 1 holds for all n G N and for ah A; > 1, see KT\i . 

K(\) K(\) 

^NxiL^^J = E NxiL^^J < J2 NxiL^,J < Y.^x{L'^,n) + 0{^). (4.7) 

fc>l k>l k>l k>l 

The latter implies ()4.5p since N\{B^) grows linearly in A when \Q\ < oo as mentioned above. 
To prove (14. 6p we consider the operators obtained from B^ by putting additional 

Dirichlet boundary condition at {x = m}, m> n. From |SV1 Thm. 1.2.1] we get 



1 1"^ 

liminf \-^Nx{B^) > liminf X~^Nx{B^J = — / fix) dx V?n > n, 

A— >oo A— s-oo ' 27r _/„ 



which implies, by letting m — )• oo, that liminf A^oo X~^Nx{B^) = oo. In view of (14.71) we 
obtain Nx [B^^ ~ Nx [B^) ■ Finally, from the Dirichlet-Neumann bracketing we deduce 
\hixtNx{B)^Nx{B^). □ 

5. Proofs of the main results 

As mentioned in the introduction, the idea of the proof is to split Nx{A(j) into two parts 
corresponding to the contribution from a finite part of and from the tail. 
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5.1. Step 1. Here we show that the contribution from the part of where x < n obeys the 
Weyl asymptotic irrespectively of the boundary conditions. Let us define 

Qn '■= {{x,y) e U : 1 < X < n} , En '■= ^ \ ^n- 

We denote by ; and Qj^,. the quadratic forms defined by the reduction of Qa on 0„ and 
En and acting on the functions from C^(r2„) and CQ{En) respectively. Moreover, let 
and Sn be the operators associated with the closures of the forms Q^i and Q^^- 

Similarly we denote by Q^^ and Qnr respective quadratic forms which are defined 
in the same way as Q^^ and Qnr but with the additional Dirichlet boundary condition at 
{x = n}. We then denote by and the operators associated with the closures of 
the forms Q^^ and Qnr- From the Dirichlet-Neumann bracketing we obtain the operator 
inequality 

r„^e5^<A^<ri^e5f, nen, (5.1) 

which implies that 

Nx{T^) + Nx{S^) < Nx{A^) < iVA(r„^) + Nx{S^), ?i G N, A > 0. (5.2) 
Lemma 5.1. For any n G'N it holds 

hm X-^NxiT^)= lim NxiTj;") =— / fix)dx. (5.3) 

A— s-oo A— s-oo zvr Ji 

Proof. Fix n G N. Since a is bounded and {— f (x) , f {x)) is for every x G (1,?^) con- 
tinuously embedded into L°°(— /(x), /(x)), it follows that there exists a constant Cn such 
that 

1 1 2 1 1 2 TV r 1 1 1 1 1 2 ( \\ 1 1 2 1 1 2 

holds for all u G C^(On)- Hence the domain of the closure of the quadratic form Q^^ is a 
subset oiH^{9.n)- The same reasoning shows that the domain of the closure of Qni contains 
the space HQ{yt,n)- From the fact that a > and from the variational principle we thus 
conclude that 

iV,(-A^J < Nx{T^) < Nx{t:^) < Nxi-A^J, (5.4) 

where — and — denote the Dirichlet and Neumann Laplacian on Qn respectively. 
Since Qn has the //^—extension property, the Weyl formula 



holds for both -A^^ and -A^^, see [Me llBiSo] . [Ng] . In view of ([^ . this completes the 
proof. □ 

5.2. Step 2. Next we will treat the contribution to the counting function N\{Aa) from the 
tail of 0. Our first aim is to transform the boundary term in (|2.ip into en effective additional 
potential. To this end we use a ground state representation for the test functions ip. Let 
/i(x) be the first eigenvalue of the one-dimensional problem 

-dyv{x,y) = fi{x)v{x,y), (5.5) 
dyv{x, -f{x)) = a{x) v{x, -f{x)), dyv{x, f{x)) = -a{x) v{x, f{x)) 
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with the corresponding eigenfunction v. By lemma lA.21 < v < 1 and v{x, y) — )• 1 as x — )• oo 
uniformly in y. Moreover, v G C'^{En). Thus every function -0 G D{Q^) can be written as 

'il){x,y) =v{x,y)ip{x,y), ip <^ Cl{En) ■ (5.6) 

Similarly, for every function %(: G D{Q^) we have 

ilj{x,y) = v{x,y)ip{x,y), ip e CQ{En) D {ip : (p{n, ■) = 0} (5.7) 

In view of (15. 6p and (j5.7p we can thus identify [V'] and Q,^ ["0] with quadratic forms [tp] 
and [ip] given by 

and acting in the weighted space L'^{En,v^dxdy). The forms Q^[(^] and Q^[ip] are defined on 
D{Qn) = C^(En) and £'(Q^) = C^(^n) n {99 : ip{n, •) = 0} respectively. A straightforward 
calculation based on integration by parts in y then gives 

Qn '""M = / ildxMf + ^^ix) \ipf + \dy^\^) dxdy. (5.8) 

J En 

Since t; — )■ 1 and /i(x) ~ a{x)/ f{x) as x ^ 00, see appendix, it is natural to compare Qn'^ 
with the quadratic form 

= / {\dM' + \dy^? + Tr\ 1^1') ^^^y- 

J En J\^) 

Let &^ and be the operators in L^{En) generated by the closures of the quadratic form 
qn[u\ on D{Q^) and D{Ql^) respectively. 

Lemma 5.2. Suppose that assumptions IKM and \2.5\ are satisfied. For any £ there exists an 
Ns such that for all n > Ns 

NxiS^)<Nx{{l-e)&^ -e), Nx{S^) > N^Hl + e)6^ + e) A > 0. (5.9) 

Proof. Let e > and let 93 belong to the domain of the quadratic forms {Qn)- From the 
fact that 

Uiix] fix] 

lim v{x,y) = 1, lim dxv{x,y) = (uniformly in y), lim — — — = 1, 

x—>-oo x— >oo X'— >oo O'ix) 

see Lemma lA. 21 and from the estimate 

\2vdxV ip dx(p\ < e\dxip\'^v^ + \ip\'^\dxv\'^ 
we conclude that for n large enough 



M > (1 - 6) 1^ (^\dM' + \dy^\' + ^ I^P) dxdy - e M^^^^ 

Q^[^] < (1 + {\dM' + \dy^\' + ^ IV'P) dxdy + e ||<^||i. 

Moreover, by ()A.6P we also have |f | < 1 so that (still for n large enough) 



) 

(En) 



(5.10) 



(1 - e)llvllL2(£;„) ^ / Ifl ^ dxdy < \ML2(^En)- 

J En 

Equation (15. 9p then follows from the variational principle by choosing e in appropriate way 
(depending on e). □ 
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5.3. Step 3. We transform the problem of studying the Laplace operator on En to the 
problem of studying a modified operator on the simpler domain 

Dn = (n,oo) X (-1, 1). 

To this end we introduce the transformation U : L'^{En) — )• L'^{Dn) defined by 



{U^){x, t) = ,/fix) ip{x, fix) t), {x, t) G Dn. 

Let An and Al^ be the operators associated with the closure of the form 

Qn[u] := qn[U~^u], (5.11) 

on CgiDn) and Cg(D.„)n{ti : n(n, •) = 0} respectively. Since U maps L^{En) unitarily onto 
L'^{Dn) and U C^iEn) = C^iDn), the variational principle gives 

N,iA:i) = N,{e^), N,{A^) = N,ie^). (5.12) 



By a direct calculation 

^\OxU\- ^ vv a u- — ' ' ' 

Now Let r/ G (0, 1) be arbitrary. Since |t| < 1 we get 



Qn[u]= / {\dM'^ + W,U^-2t^-d:,udtU+^—{tudtU + t^\dtu\'^) + ^\dtu\'^)dxdt. 

J Dr, J J J 



— dxudtu\ < r] Idxul"^ + 1]^^ Idtul"^ , -j^\tudtu\ < r]^^ \u.\ -r 



2 , '/ lfi,.|2 



Moreover, from (|2.4p and from the fact / is decreasing at infinity, by assumption 12.21 it 
follows that 

f'{xf<2f{x) sup \ns)\, (5.13) 

s>x 

for all X large enough. Since /" — )• as x — )• oo, for any tj £ (0, 1) there clearly exists an A'^ 
such that for any n > N^j it holds 

Qn[u]^ [ {il±v)\dM^ + W„u^ + ^^^\dtu\^±r]u^)dxdt. (5.14) 

JD„ J 

We denote by Hn and Hn the operators acting in L'^{Qn,r) associated with the closures of 
the quadratic form 

/ {\d'M^ + Srr + ^-(^) dxdt 

JDr^ f [X] 

defined on C'^{Dn) and CQ{Dn) H {u : u{n, •) = 0} respectively. 

Lemma 5.3. Suppose that assumptions \2.2\ and \2.5\ are satisfied. For any e there exists an 
such that for all n > and any X > it holds 

N,{A^)<N,iil-e)H:^), N,{A^)>N,i{l + e)H^). (5.15) 

Proof. In view of the fact that Wa{x) — t- oo the statement follows from (j5.14p . □ 

Next we observe that since depends only on x, the matrix representations of the operators 
Hn and Hj^ in the basis of (normalised) eigenfunctions of the operator —f{x)~^ ^ on the 
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interval (—1, 1) with Neumann boundary conditions are diagonal. We thus have the following 
unitary equivalence: 

oo oo ,2 1.2 2 

i^n^-e<n, ^f-e^?., <'' = -|2 +^'^(-) + ^' (5-16) 

where H^!^ are operators in L'^{n, oo) with Neumann/Dirichlet boundary condition at x = 
n. We denote 

njN _ njN njD _ njD 
'T-O.n — i^n 1 "-O.n " "-n • 



As a consequence of (j5.16p we get 

Proof of Theorem 13. it We make use of inequality (jS.ip for some fixed n and show that the 
operator on the left hand side of ()5.ip has purely discrete spectrum. By general arguments 
of the spectral theory this will imply the statement. Since the spectrum of is discrete, 
it suffices to show that the same is true for . In view of Lemma [5. 21 and equations ()5.12p . 
(|5.15p it is enough to prove the discreteness of the spectrum of . By (|5.16p we have 

spect(/7^) = U^=ospect(1^^J, 

First we notice that spect('H^^) is purely discrete for each k and n. Indeed, a sufficient 
condition for the spectrum of to be purely discrete is that 



^o-(a;) + 77732 — ^ C)0 as x — )■ oo , (5-17) 



4/(x)2 

see e.g. |RSl Thm. 13.67], which is a direct consequence of assumption (|2.5p . Hence the spec- 
trum of is pure point, i.e. consists only of eigenvalues. Moreover, since f'^{x)Wtj{x) — )• 
as X — )• oo by ()5.13p and boundedness of o", it is easy to see that 

Vn inf spect(?^^„) — )■ oo as /c — )• oo. 

Hence all the eigenvalues in the spectrum of have finite multiplicity and spect(i?^) 
contains no finite point of accumulation. This means that spect{H^) is discrete. □ 

Proof of Theorem \3.3[ Case < oo. By Lemma [4. li the asymptotic behaviour of N\{T-Ln'^) 
does not depend on the boundary condition at x = n, nor on n itself: 

Nx{n^) ^ Nx{n^) ^ Nx(n^) A^oo, VneN. (5.18) 

Now fix an e > 0. From Lemma 15.21 (j5.12p and (|5.15p we see that for all n large enough it 
holds 

Nx{S^) < Nxiil - e) H^). Nx{S^) > iV,((l + e) H^) (5.19) 
Moreover, f'^{x)Wa-{x) — t- at infinity so that 

, . k TT k 7T , , k 7T , 

(l-^>47(iP^'^'(^) + 47W^S(l+e)j^ VA->1 (5.20) 

for all X large enough uniformly in k. Now observe that the sequence {A;27r2/4/(x)2}^^^ 
enlists all the eigenvalues of the operator — /(x)"^ ^ on the interval (—1,1) with Dirichlet 
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boundary conditions. Hence it follows from ()5.16p and (I5.20p that for n large enough 
iV,((l + e) H^) > N,{{1 + ef B^) + iV,((l + e) U^) 
Nx{{l-e)H!!)< N,{{l-efB:!) + N,{{l-e)n^), (5.21) 

where Bn^^ are the operators defined in Section HI Note that Bn^^ and Hn^^ satisfy 
assumption (I2.10p by Lemma 14.51 and equation ()5.18p . In view of ()5.2p we then conclude 
that for n large enough 

NM.) < (1 + 0{e)) {N^{T^) + N,{B^) + N,{n^)) (5.22) 
N,{A^) > (1 + 0{e)) {N,{T^) + N,{B^) + N,{n^)) , (5.23) 

If the volume of is finite then it follows from Lemmas 14.51 15.11 and equations (j5.18p , (j5.22p , 
(fOajl that for any e > 

1 + 0(e) < liminf , ^^(^a) ^ j.^ ^ Nx{A,) ^ ^ ^ ^ 

A|0|+iV;,(7^,) - A^oo ^|l^|+iVA(^.) - 



By letting e — )■ we arrive at (|3.ip . 

Case |fi| = oo. If the volume of il. is infinite, then Lemma 14.51 gives 

Nx{T^) + Nx{B^) ~ Nx{B^) ~ Nx{B^) ~ iVA(Tf ) + NxiB^) ~ iVA(^) (5.24) 

as A — 7- oo. Moreover, mimicking all the above estimates for the Dirichlet-Laplacian — 
instead of A^^ it is straightforward to verify that for any e > and n large enough, depending 
on e, it holds 

Nx{{l-e)B^) < Nx{-AE) < Nx{il + e)BS). 

This together with (j2J0|) and (jOi]) implies that A''a(^) ~ Nx{-A^) as A oo. Equation 
dHU thus follows again from K22li and (fOHjl . □ 



Proof of Corollary \3.8[ Note that the assumption 12.51 is fulfilled. Indeed, equation ()5.13p 
shows that f{x)V{x) — )• 0. Consequently (j2.6p holds true since / — )■ and 



W^{x) ^ ^ X ^ oo. (5.25) 
fix) 



To prove ()3.3p we recall that 



liminf A-^iVA(A^) > liminf X'^Nx{-A^) = 



On the other hand, if limsup^^o^ x^/(x) = 0, then (|5.25p says for any e > exists an Xg 
such that W(j{x) > p- holds for all x > x^. Lemma 14.41 gives 

-l2 ^2' 



/ (f x^\ 
limsup X'^NxCHa) < limsup A^^iVA -y^ + — 

A— 5>oo A— >cxD \ ax £ J j^2u 



e 
4 



From Lemma 14.51 and the proof of Theorem 13.31 see equations (|5.2p , (j5.3p , ()5.19p and (j5.2ip 
we then get 

limsup ^ — - < (1 + 0[e)) + limsup — - — < (1 + 0{e)) + 0[e) . 

A^oo A 47r A^oo A 47r 
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Equation (13. 3p now follows by letting e — )• 0. In order to prove (I3.4p we note that Wa{x) ~ 
a a'"^ as X — >■ oo, see ()5.25p . From Lemma [13] we thus deduce that 

lim X-^NxCHc 



SO that (|2.9p is satisfied and (I3.4p follows from Theorem 13. 61 □ 

6. Generalisations 

6.1. Non symmetric boundary conditions. As mentioned in Remark 13.51 the above 
approach can be applied also to Robin Laplacians with different boundary conditions on the 
upper and lower boundary of Q. More precisely, to operators A^-^^^^ generated by the closure 
of the form 

Qai,a2[u] = / \Vu\'^dxdy+ / {ai{x)u{x,f{x))'^+a2{x)u{x,-f{x)f)dx (6.1) 
Jn Ji 

on Cq{Q). We can proceed in the same way as in section [5] replacing the function v{x,y) in 
step 2 by the function w{x,y), which solves the eigenvalue problem 

-dyw{x,y) = ll{x)w{x,y), (6.2) 
dyw{x, -f{x)) = cri(x) w{x, -f{x)), dywix, f{x)) = -a2{x) w{x, f{x)), 

with p.{x) being the principle eigenvalue. From equation (jA.lip . see appendix, we then get 
a generalisation of Theorem 13. 3i 



Proposition 6.1. Suppose that assumptions \2.5\ and \2.6\ for ui . o"2 are satisfied. Then 

Nx{A^,^.,) ~ A^a(-A^) + Nx{n,) A ^ oo, a{x) = ^^i'^) + (5.3) 

6.2. Dirichlet-Neumann Laplacian. Our second remark concerns the case in which we 
impose Dirichlet boundary condition on one of the boundaries of 0. We confine ourselves 
to the special situation when we have Dirichlet boundary condition on one boundary and 
Neumann on the other. We denote the resulting operator by ^0,00 • 

Proposition 6.2. Let < 00 and assume that f is decreasing at infinity. Then 

NxiAo^oc) ~ P A, A ^ 00. (6.4) 
47r 

Proof. First we observe that by the variational principle. 

liminf A-^iVA(^Ooo) > hminf X'^Nx{-A^) = . (6.5) 

A— >oo ' A— >-oo AtT 

Assume that / is decreasing on (a, 00) and that A is large enough so that there exists a 
unique point x\ > a such that f{x\) = 7r/(4-v/A). We impose additional Neumann boundary 
condition at {x = x\} dividing thus O into the finite part Oa := {{x,y) £ : x < xx} 
and its complement il^. It is then easy to see that the quadratic form of the corresponding 
operator acting on Qj^ is bounded from below by 

/ ( , „ .0. X + \dxu\'^] dydx > A / / u dy dx 
XX J -fix) V16/^(x) J J^f^^x) 
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for all functions u from its domain. Consequently, this operator does not have any eigenvalues 
below A. To estimate the number of eigenvalues of the operator acting on ilx^ we cover 
with a finite collection of disjoint cubes of size L = l/{eVX) with e > 0. Since ^Ix has the 
extension property, the standard technique of Neumann bracketing gives 

X^'NxiAo,^) < X-'Nxi-Al) < ^ (1 + 0{s)) + c ^ (1 + e"^) 

< A-iiV,(-A^J< ^(i + 0(e)) + c^(i + e-i), (6.6) 

where c is independent of A. However, since / is integrable and decreasing at infinity it is 
easily seen that x/(x) — )• as x — oo. Hence 

x\ 4 

limsup —= = — limsup xxf{x\) = 0. 

A— >oo V A ^ A— >oo 

Letting first A — )• oo and then e — )• in (j6.6|) we obtain limsup;^„^o^ A~^A^a(^o,oo) < l^l/^vr, 
which together with (j6.5p implies the statement. □ 

Appendix A. 

Lemma A.l. Let n{x) be the function defined by the problem (jS.Sp . Then 

Mx) < ^ Vx>L (A.l) 
fix) 

Proof. For each fixed x € (1, oo) we define the quadratic form 

aM= \u'{y)\^dy + a{x){\u{f{x))\^ + \u{-f{x))\^), u € D{a,), (A.2) 
where D{ax) = {— f (x) , f (x)) . The variational definition of /i says that 
M(-) = mf ^ „_ < - ^ 



□ 



In the next Lemma we use the notation k{x) := ^ l^ix). 



Lemma A.2. Let the assumption \2.5\ be satisfied. Then the eigenfunction v{x,y) of the 

problem (15. 5p associated to the eigenvalue f^{x) is twice continuously differentiable in x. 
Moreover we have 

lim = 1 (A.3) 

lim v{x, y) = 1 uniformly in y, (A. 4) 

lim dxv{x, y) = uniformly in y. (A. 5) 

Proof. It is easy to see that 

v{x,y) = cos{K{x)y), (A.6) 
where k{x) is the first positive solution to the implicit equation 

F{x, k) := K tan(«;/(x)) - a{x) = 0. (A.7) 
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Since f{x)K{x) — )• as x — )• oo by Lemma [A. II (recalling that a{x)f{x) — )• 0), we easily 
deduce from ()A.7p that 

f{x)ti'^{x) ^ ^ 



lim ' '-^f^-^' = 1, (A. 



which proves ()A.3p . Equation ()A.4p thus follows directly from ()A.6p and the fact that 
/(x)k(x) — )• 0. Next we note that ()A.7p implies 



and hence 



d^F{x, = tan(/(x)A^) + {^''J \ > 0. (A. 

cos [J[X)K) 



Since a S C^(l,oo), the implicit function theorem shows that k is of the class and in 
view of ()A.6P we see that f is twice continuously differentiable in x. 

In order to prove (jA.Sp we need some information about the behaviour of k' for large 
From the positivity of / and a and from the Taylor theorem we conclude that a' /y/a is 
bounded and that f / y/f — )■ 0, see equation (j5.13p . Equations (\A.9\i and (jA.SP then give 



k'(x) = ~ (/'(x) - x^oo. (A.IO) 



On the other hand, a direct calculation shows that 



\dxv{x,y)\ < \k'{x)\ f''^{x)^/^ Vx>l. 
This implies (IXSl) . □ 

Notice that if we replace the eigenvalue problem (|5.5p by (j6.2p . then a straightforward 
analysis of the associated implicit equation shows that 

hm_ — ^ = 1, a{x) = . (A.ll) 
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